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ADIABATICS  OP  SUPERHEATED  STEAM 


I  INTRODUCTION 

The  rapid  development  in  the  use  of  superheated  steam 
within  recent  years  has  made  it  highly  desirable  that  there 
should  be  available  an  easy  means  for  making  calculations 
which  arise  in  engine  and  turbine  practice.    All  of  the  form- 
ulae now  available  for  making  these  calculations  are  exceed- 
ingly difficult  to  apply  in  practice,  and  further,  the  con- 
flicting nature  of  the  of  the  experimental  data  upon  which 
they  are  based  cast  doubt  upon  their  reliability*    It  would 
seem  that  tables  of  the  properties  of  superheated  steam 
such  as  those  by  Professor  Peabody  and  by  Professor  Marks 
and  Dr.  Davis  would  simplify  such  calculations,  but  these 
again  are  subject  to  the  objection  that  they  were  obtained 
by  the  graphical  integration  of  assumed  specific  heat  curves. 

In  a  recent  paper  by  Professor  Goodenough  a  connected 
theory  has  been  evolved,  by  means  of  which  the  important 
properties  of  superheated  steam  can  be  calculated  directly 
without  recourse  to  approximation.    This  theory  is  based  upon 
the  most  reliable  experimental  data  available. 

(a)      EQUATIONS  FOR  SUPERHEATED  STEAM. 

The  following  derivation  of  equations  is  taken  from 
the  above  mentioned  paper  by  Professor  Goodenough. 

"Prom  the  general  principles  of  thermodynamics  we  have 
the  relation 


2 


6  Cj 
6P 


=  -AT 


St2 


(i) 


which  applies  to  all  substances.  Writing  the  characteristic 
equation 


pv  =  BT  - 


(  x  +  ap  )-£-  p, 


in  the  form 


(2) 


(3) 


we  readily  obtain  "by  differentiation, 


&  v  _  B      (  -i      „_  }  mn 


■f^=-(  1  +  ap  )  S3  <n  + 
£  ipn  +  -s 


whence  for  superheated  steam  (1)  becomes 


6  «, 


=  (  1  ♦  ap)  ^  (n  *ll 

<j01  ♦  1 


(4) 


Taking    T    as  a  constant  and  integrating  (4)  with    p    as  the 
independent  variable  the  result  is 

cp  =  (  1  ♦  |p  )p  ^SSlJM^JsI  +  constant  of  integration. 

The  constant  of  integration  may  be  a  function  of    T,    since  T 
was  held  constant.  Hence 


c    =  <D(T)  +  (  1  *  *p  )  (n  »  1) 

P  2  rpn  ♦  1 


(5) 


or 


C    =  <D(T)  +  f  (p,T),  --(6) 

P 


if  the  term 


,  _  a  „\  Amn  (n  +  1)  w 
<  1  +  2p)       ^  ♦  1  p 


is  replaced  for  convenience  by  f(p,T)S 

Prom  a  consideration  of  the  most  reliable  experimental 
data  available,  Professor  Goodenough  shows  that  the  function 
(D(T)    has  the  form  <*  for  temperatures  usually  employed 

in  engine  and  turbine  practice.    Whence  ,  substituting  in  (5) 

Cp=  (  i  *  |  p) js^hnr*  » <7> 

"With  the  aid  of  an  explicit  equation  for  the  specific 
heat,  expressions  for  the  heat  content,  energy  and  entropy  are 
readily  derived.    For  this  purpose  the  general  equation 

dq  =  cp  dT  -  AT  |^      dp      -  (3) 

is  most  convenient.    Since  heat  content  is  defined  by  the 
relation 

i  =  A  (  u  ♦  pv  ) 

we  have 

di  ■  dq  ■#■  Av  dp, 
whence  from  (8) 


di  ■  e    dT  -  A(  T  &  -  v)  dp 
P  6t 


(9) 


t 
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From  the  characteristic  equation  (3)  the  derivative  is 

<5  T 

found  to  be 

*&  v  _    B  v  mn 

-~  =    =•  *  (  1  ♦  ap  )  •  • 

Substituting  this  and  also  the  expression  for  cp  in  (9),  the 
result  is 

di  =  (*  dT  4  Amn  (n  ♦  1)  p  (1  +    §  P^-y  + 

Amn  (n  ♦  1 )_  /,  ,^ 
 ■  — ■*  CI  ♦  ap  )  dp 

The  integral  of  this  exact  differential  is 


i  =    T+iTS-  fa.  t  li  p(i  +  &  p)  »  i   (10) 

2 


"Since    Au  =  i  -  Apv,  we  obtain  for  the  energy  expressed 
in  thermal  units 

Au  *  T(*  %    f   T2-  AB)  -    ~L[n  ♦  (n  +  1)  |  p]  +  iQ  (11) 


"For  the  entropy    s    an  expression    is  readily  found 
from    (8)«  Thus 

dS  =       =  CPT  "  A  S  dp   (12) 


* 


Introducing  in  this  equation  the  expression  previously  derived 
for    Cp    and     ^i,  the  result  is 

dB  s  (|  ♦    /?)dT  +  Amn(n  +  1)  p  (1  +  §  p)^d»  g    -    AB  S|  . 

(1  +  ap)  dp  (13) 


T11 


The  integral  is 

s  =   °<logeT  +x?T  -  AB  logep  -  1(  1  ♦    |  p)  ♦  b0    -  -  (14) 

The    constants  in  the  above  equations  have  the  following 
values  as  determined  by  Professor  Goodenough. 

English  units* 

gL  «  0.372 

jS  =  o.oooi 

log  m  ■  9.50308 
n  *  3.5 

4-    — I 
777.7  J 

T  =  t  +  459.6 

B  =  85.77     (  p  in  pounds  per  square  foot) 

*  0.5956  (  p  in  pounds  per  square  inch) 

iQ  =  882.4 

s     =  -0.4300 
o 

a  =  0.0007 
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Reducing  to  common  logarithms,  and  suds titu ting  the 
above  constants,  equation    (14)    reduces  to  the  form 

8  ■  0.85657logT  +  0.0001T  -  0.253921ogp  -    pff*  1  +  |  p)  - 

0  *  4300   (15) 


(b )     ADIABATIC  EXPANSION  OF  SUPERHEATED  STEAM. 

A  problem  frequently  met  with  in  steam  engineering  is 
that  of  calculating  the  work  done  during  the  adiabatic  expan- 
sion of  steam,  since  the  adiabatic  expansion  represents  most 
nearly  the  actual  change  taking  place  in  the  cylinder  of  an 
engine  or  the  stages  of  a  turbine.    It  is  highly  desirable, 
therefore,  that  a  simple  method  of  making  these  calculations 
should  be  provided. 

By  the  use  of  equation    (14)  we  are  enabled  to  find  the 
pressure  corresponding  to  any  temperature  during  an  adiabatic 
expansion,  since    s    is  then  constant,  and  by  substituting 
corresponding  values  of    p    and    T    in  the  modified  form  of 
the  characteristic  equation,     (3),  the  corresponding  value  of 
v,     the  volume  may  be  determined.    However  this  is  an  exceed- 
ingly laborious  process,  involving  as  it  does  the  solution  of 
the  transcendental  equation    (14).    It  is  evident  that  if  an 
exponential  relation  between  pressure  and  volume  of  superheat- 
ed steam    could  be  found  similar  to  that  for  a  perfect  gas, 
namely      pv11  =  C,    and  to  that  deduced  by    Zeuner    for  satur- 


7 

ated  steam,  the  labor  of  calculating  volumes,  and  work  done 
during  expansion  would  "be  materially  decreased. 

It  is  the  purpose  of  the  writer  to  derive  if  possible 
an  exponent  that  will  give  results  with  an  accuracy  consistent 
with  that  of  the  data  generally  available  in  engineering 
problems. 
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II      METHOD  OP  DETERMINING  THE  EXPONENT. 

Equation    (15)    will  "be  taken  as  the  "basis  of  this 
investigation.    Selecting  some  convenient  value  of    s    on  the 
T-s    chart,  (for  this  purpose  the  chart  accompanying  the  Marks 
and  Davis  Steam  Tables  is  convenient)  the  chosen  value  of 
entropy  is  substituted  for    s    in  equation    (15),  which  then 
reduces  to  a  transcendental  with  the  two  variables    p    and  T. 
If  either    p    or    T    is  given  an  assumed  value  the  other  may 
be  calculated.    It  seems  most  convenient  for  our  present 
purpose  to  assume  values  of    T    and  calculate  the  correspond- 
ing value  of    p.    For  a  given  value  of    s    and    T,  equation 
(15)  may  be  written 

logp  M-K  (pi    |  P2)    -  -   -(16) 

where  (0.856571ogT  t  0.0001T  -  s  -  0,4300 )/0. 25392  ■  C, 
and 

m  n  A        n  K< 

0.25392T11  +  1 

The  solution  of  equation  (16)  is  accomplished  most  convenient- 
ly   by  writing  it  in  the  form 

logp  =  y  =  C-  K(p  +    |  p2), 

and  plotting  the  two  curves  as  shown  below. 


r 

Evidently  the  intersection  of  these  curves  gives  the  value  of 
p    that  satisfies  equation  (16).    As  a  check  upon  the  solution, 
the  value  of    p    thus  obtained  is  substituted  back  in  equation 
(15).    In  a  similar  manner  other  values  of    p    are  obtained  on 
the  same  adiabatic  by  using  different  values  of    T.    In  these 
calculations  the  value  of    T    varies  from    750°  to    1500°  by 
increments  of    100°  .    The  lower  points  being  taken  near  the 
saturation  line,  and  the  upper  points  are  higher  than  temper- 
atures employed  in  practice. 

Having  values  of  p  and  T  on  an  adiabatic,  the  corre- 
sponding values  of  v  may  be  found  by  substituting  p  and  T 
in  equation  (3).  Having  a  series  of  pressures  and  volumes  on 
an  adiabatic  it  becomes  possible  to  determine  the  relation  be- 
tween them. 

This  same  process  was  continued  for  eight  different 
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values  of  entropy,  varying  from    1.60    to    1.918,  covering  the 
regions  usually  employed  in  practice. 

Having  pressures  and  volumes  for  different  entropies,  it 
becomes  possible  to  determine  whether  or  not  the  relation  of  p 
to    v    is  the  same  throughout  the  regions  under  consideration. 

To  determine  the  exponent    n    in  the  expression    pv11  =  C, 
which  it  is  intended  to  establish,  pressures  and  volumes  for 
each  entropy  used  are  plotted  on  logarithmic  cross -section 
paper  as  shown  in    Figure  1.    The  slope  of  the  line  drawn 
through  the  points  corresponding  to  a  particular  entropy, gives 
the  desired  exponent.    The  slope  is  determined  by  scaling  the 
intercepts  of  the  lines  with  the  axes. 
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III       RESULTS  MP  C ONCLUSI ONS . 

Eight  values  of  entropy  were  used  in  these  calculations, 
and  the  value  of  the  exponent    n    corresponding  to  each  is 
given  in  table  I. 


TABLE 


Entropy 
1.6000 
1.6571 
1.6726 
1.70459 
1.7248 
1.74809 
1.80507 
1.91800 
Mean  = 


Value  of  n 

1.301 

1.303 

1.300 

1.293 

1.294 

1.298 

1.283 

1.290 

1.295 


Then  we  may  write 


1.295 


also 


pv 


v2  =  vx 


=  C, 


Pi 
l?2J 


1  

1.295 


and 


W  =    plvl"  p2v2 


1.295  -  1 


x  144 


(if    p    is  in  pounds  per  square  inch)    for  the  work  in  foot 
pounds  which  is  done  during  adiabatic  expansion  from  a  pressure 


12 


p^  to  a  pressure  pg. 

As  a  check  upon  the  accuracy  of  this  exponent,  the  work 
in  foot  pounds  is  in  each  case  compared  with  the  total  change 
of  energy  as  given  by  the  expression 


U2  =  T2(  202***  *  0.0389T2)  -  +  0.00025  p8)  +  686242 


Log  c  =  12.20551    (  p    in  pounds  per  square  inch  ) 
If  the  exponent  just  derived  for  the  adiabatic  expansion  has 
the  proper  value,    W  =       -  U2,    where       is  the  total  energy 
ifl  foot  pounds  at  the  higher  temperature  and  pressure  and 
Ug  is  the  total  energy  in  foot  pounds  at  the  lower  temperature 
and  pressure.      The  error  in  per  cent  of    W    as  compared  with 
U    -  U      is  given  in  table  II.  In  table    II    are  given  also 
the  upper    and  the  lower  pressures  and  volumes  used  in  mak- 
ing the  calculations;  also  the  lower  volumes  as  calculated  by 
using    the  relation    P^v^  =  Pgv§,    Pi  »  P2>  811(1  vl  fcQing 
known.    The  error  in  per  cent  of  this  approximate  value  of  v2 
as  compared  with    Vg  calculated  by  means  of  equation    (3)  is 
also  given. 


Ux  m  Tx(  202.44  +  0.0389^) 


CPl 


(1  +  0.00025  p.,  )  +  686242 
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It  will  "be  observed  from  table    II    that  the  value  1.3 
for  the  exponent  gives  results  all  of  which  are  accurate  to 
within  less  than  one  per  cent  when  used  at  the  entropies  1.65 
to  1.75,  which  is  the  region  usually  covered  in  actual 
practice.    At  the  higher  entropies,     1.8    to    1.9  ,  the  agree- 
ment is  not  so  close  when  using  this  value  for  the  exponent. 
At  these  higher  entropies  the  value  1.295  gives  more  accurate 
results. 

This  exponent  ought  in  all  probability  to  be  some  func- 
tion of  the  entropy  as  will  be  observed  from  table  I.  How- 
ever this  would  again  complicate  matters,  and  thus  defeat  the 
purpose  of  this  paper.  It  seems  advisable  therefore,  to  take 
the  value  1.3  as  representing  nearly  enough  for  all  practi- 
cal purposes  the  exponent  of  the  adiabatic  curve  of  super- 
heated steam. 
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